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AN  IMPROVED  SCALING  LAW  FOR  DETERMINING  STIFFNESS  FOR  FUT. 

CIRCULAR  CANOPIES 


INTRODUCTION 

In  the  design  and  developnient  of  parachutes,  it  is  often  convenient  to  rely  on  scale 
models  for  drop  testing  or  wind  tunnel  testing.  The  use  of  scale  models  in  development  testing 
can  save  significant  cost  and  time  in  the  development  cycle.  Because  the  parachute  scaling 
laws  are  not  completely  understood,  the  interpretation  of  model  test  data  posss  numerous 
problems.  For  example,  model  parachutes  tend  to  open  quicker,  with  a  larger  scaled  opening 
shock  force,  than  their  full-scale  counterparts. 

For  the  scaling  of  parachutes,  there  are  several  unsettled  questions  concerning  the 
scaling  of  the  fabric  stiffness.  Scale  models  of  parachutes  have  customarily  been  made  from 
the  same  ^  iight  fabric  as  the  prototype.  This  makes  the  scale  model  canopy  relatively  stiff 
compared  to  full-scale.  Heinrich  recognized  this  and  developed  a  stiffness  index  to  compare 
model  and  iimtotype  and  did  research  on  canopy  construction  techniques  needed  to  make  the  model 
fabric  stiffness  more  nearly  that  of  the  prototype.  However,  even  using  modified  construction 
techniques,  it  has  not  been  possible  to  find  ways  to  make  the  models  quite  as  flexible  as  the 
prototype. 

This  report  presents  an  investigation  of  the  parachute  scaling  laws  applicable  to  the 
opening  process.  The  object  w’as  to  find  a  better  correlation  parameter  for  opening  time  so 
models  can  still  be  suffer  than  the  prototypes,  but  that  model  test  conditions  can  be  altered 
to  make  the  model  data  more  usable  to  predict  full-scale  performance.  To  accomplish  this, 
existing  theories  for  parachute  opening  dynamics  were  reviewed,  existing  test  data  were 
analyzed,  and  parachute  canopy  fabric  properties  were  examined.  No  new  theories  were  developed 
and  no  new  experiments  were  conducted  as  part  of  this  study.  The  end  result  of  this  study  was 
the  identification  of  a  better  correlation  parameter  for  opening  time  than  any  currently 
available  to  predict  prototype  parachute  performance. 


REVIEW  OF  DIMENSIONAL  ANALYSES  AND  SCALING  PARAMETERS 


Over  the  yetis,  t  number  of  dimensional  analyses  of  the  parachute  filling  process  have 
been  conducted.  These  analyses  have  identified  the  importance  of  the  Froude  number  (ratio  of 
inertia  to  gravity  forces)  and  mass  ratio  on  parachute  opening  time  and  shock.  Other 
parameteis  identified  as  being  of  somewhat  lesser  importance  are  the  Reynolds  number,  fabric 
porosity,  shroudline  elasticity  and  damping,  canopy  stiffness,  and  canopy  mass  ratio  (or 
inertia  number). 

The  more  important  of  these  investigations  have  been  conducted  by  Kaplun 
[1],  Barton  (2],  French  [3,4],  Fu  (5),  and  Lingard  (6).  Only  a  summary  of  the  major 
results  wnll  be  presented  here. 

One  of  the  first  investigators  to  point  out  the  importance  of  parachute  elastic 
properties  was  Kaplun  (1).  Quoting  from  his  work: 

The  specification  of  the  elastic  properties  of  the  parachute  is  so  complicated  that,  at  least, 
some  departures  from  rigor  wili  be  unavoidable.  Properly,  considering  the  fabric  to  be  a  plastic, 
anisotropic  medium,  it  would  be  necessary  to  specify  twenty>one  elastic  parameters,  and,  in 
addition  to  the  similarity  of  distribution,  require  that  the  stress>strain  diagrams  remain 
similar.  Let  us  tentatively  take  the  other  extreme  by  assuming  that  the  spring  constant  of  the 
shroudlines  and  the  flexural  rigidity  of  the  canopy  are  the  only  important  factors.  Whether 
this  assumption  is  actually  justified  can  be  decided  only  on  the  basis  of  consistency  of  the 
experimental  data  derived  from  it.  However,  it  should  be  mentioned  that  rigidity-,  which 
vanes  inversely  as  the  fourth  power  of  the  diameter,  is  especially  important  in  small  models. 


Kaplun  considered  the  case  of  a  parachute  attached  to  a  fixed  support,  thus  the 
Froude  number  did  not  appear  in  the  results  of  his  evaluation.  He  came  up  with  a 
set  of  six  dimensionless  parameters,  choosing  the  combinations  of  parameters  based 
on  expected  physical  significance.  The  list  of  parameters  selected  was  the  following: 


Reynolds  number 

Reynolds  number  of  the  fabric 

Vdp 
Rf  * - 

V 

Mach  number 

M 

Shroudlin**  elasticity 

k 

*S  *  ..1^ 

' 

pV2D 

2 


Canopy  rigidity 

Canopy  inertia  number 

Based  on  these  independent  parameters,  the  dimensionless  dependent 
parameters  were  selected  as  dimensionless  time, 

Vt 

D 

where  total  filling  time  T  is  given  by 

VT 

(R,  Rf,  M,  es,  ec,  I) 

and  dimensionless  force 


Cc 

I- 


(ED 

'pV2D^ 

m 

pD3 


f 


where  maximum  opening  force,  Fnia*,  is  given  by 


pV2D2 


f2  (R,  Rf,  M,  es,  ec,  D 


Although  Kaplun  did  not  use  Froude  number  as  one  of  his  variables,  he 
defined  it  as  follows; 


Fr  • 


5PV2D2 


Using  this  definition  with  the  case  of  a  parachute  tied  to  a  fixed  support 
(infinite  mass  inflation)  gives  a  Froude  number  of  zero.  Other  investigators  define 
Froude  number  differently,  as  pointed  out  later. 
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Kaplun  pointed  out  that,  if  it  is  assumed  that  these  parameters  have  no  effect, 
then  the  expressions  for  the  performance  variables  would  only  be  determined  up  to 
a  constant  of  proportionality,  i.e., 

VT 

*  constant 


and 


■: - =  constant 

^pV2D2 

For  certain  ranges  of  operation,  this  t\’pe  of  behavior  has  been  observed  for  the  non- 
dimensional  filling  time,  as  noted  later. 

It  should  be  noted  that  Kaplun  did  not  define  how  the  "flexural  rigidit}’,"  El, 
of  a  parachute  canopy  should  be  interpreted  or  determined,  either  analytically  or 
experimentally.  The  canopy  flexural  rigidity  must  be  evaluated  in  some  manner 
that  is  physically  consistent  with  its  effect  on  parachute  opening  d\Tiamics.  Heinrich 
[7]  did  work  in  this  area  and  defined  a  "canopy  stiffness  index"  quite  different  from 
Kaplun's  canopy  rigidity  parameter.  However,  Heinrich's  index  is  useful  only  for 
comparing  canopies  or  in  attempting  to  construct  scale  model  canopies  with  the 
same  flexibility’  as  a  prototype  canopy.  But  it  does  not  let  someone  use  data  on  a  stiff 
canopy  model  to  predict  the  performance  of  a  more  flexible  prototype.  This 
determination  is  one  of  the  contributions  of  the  current  study. 

Fu  (51  developed  a  theory  to  predict  parachute  opening  behavior,  and 
reviewed  dimensional  analyses  conducted  by  previous  investigators.  In  addition  to 
scaling  geometric  parameters  between  model  and  prototype,  Fu  identified  the 
following  parameters  as  needing  consideration: 

Vs2 

Froude  number,  Fr  «  “pr 
gt-'o 

Mass  ratio,  m  *  — T 
pDo 

Initial  trajectory  angle,  Yo 
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ncfDo 

Elastic  factor,  = - 2 

mLV, 

^  ndDo 

Damping  factor, 


Effective  porosit)’,  Cpfi' 


Payload  to  system  mass  ratio, 

The  effect  of  canopy  stiffness  was  not  considered  directly  by  Fu,  but  some  of  its 
implicit  effects  were  accounted  for  as  an  initial  input  to  the  theory,  in  the  form  of 
semiapertural  angle,  Po  (see  Fig.  1).  This  angle  determines  the  initial  parachute 
mouth  area  available  for  air  capture.  Empirical  values  for  Po  were  obtained  by  Fu 
from  films  of  parachute  opening,  and  the  value  of  Po  would  certainly  be  affected  in 
part  by  canopy  stiffness. 


In  [8],  the  important  dimensionless  correlation  parameters  are  summarized  as 
follows: 


Name 

Parameter 

Svmbol 

Factor 

Reynolds  No. 

hpip 

Rc 

Scale 

Mach  No. 

t'/Cs 

M 

Compressibilit)’ 

Froudc  No. 

Ft 

Fluid  displacement 

Kaplun  No. 

i  K/pV^I 

Kp 

Stretching  resistance  (1 /elasticity) 

Relative  stiffness 

Ellpv^l^ 

n 

Stiffness  (1 /flexibility) 

Mass  Ratio 

pl^lrr. 

Rm 

Added  air  mass 

Effective  Porosity 

vx/v‘ 

Ceff 

Air  permeability 

H 

ere  again,  the  relative  stiffness  parameter,  Q  =  — —  ,  is  presented  without  any 
dication  of  how  to  evaluate  the  canopy  “flexural  rigidit}',"  EL 
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In  [9],  Lee  shows  that  time  to  peak  opening  shock,  x,  and  the  value  of  peak 
opening  shock,  Xp,  are  functions  of  mass  ratio,  Froude  number,  canopy  stiffness 
index,  and  parachute  spring  constant: 


Xp,T-f(Mr,Fr,ii,Kp) 


where  canopy  stiffness  index  ii 


spring  constant  Kp 


ksDo 
mig  • 


It  was  shown  that  matching  of  Mr  and  Fr  and  satisfying  geometric  similarity  would 
not  necessarily  guarantee  scaling  t  and  Xp  because  other  parachute  properties  had  to 
be  considered  also;  these  include  fabric  mass  density  and  porosity,  canopy  flexibility, 
and  suspension  line  stiffness. 

Steeves  [10]  showed  from  dimensional  analysis  that  the  scaling  of  a  parachute 
opening  requires  that  the  canopy  fabric  areal  density  ratio  between  model  and  full- 
scale  be  the  same  as  the  geometric  dimensional  ratio  and  that  the  fabric  porosity 
ratio  be  the  square  root  of  the  geometric  ratio. 

Lee  concludes  that  "...the  present  moderate  amount  of  data  show  that  Xp  and 
X  of  full-scale  parachutes  at  a  given  Fr  and  Mr  can  be  reasonably  well  predicted  from 
the  measurement  of  Xp  and  x  of  small-scale  models  with  similar  t)  and  Kp  as  those 
of  the  full  scale." 

One  of  the  difficulties  in  satisfying  the  above  requirements  is  to  be  able  to 
build  small-scale  models  with  sufficiently  low  stiffness  indices,  t|,  to  match  the 
prototvpe  parachute  stiffness.  A  way  to  avoid  the  requirement  for  matching  t| 
would  be  to  develop  a  "relative"  or  "effective"  stiffness  index  that  would  allow 
scaling  data  from  stiffer  models  to  more  flexible  prototypes  by  use  of  a  suitably  ntodi" 
fied  scaling  parameter. 
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REVIEW  OF  OPENING  DYNAMICS  THEORIES 


It  was  not  the  puTpose  of  this  study  to  develop  any  new  theories  for  parachute 
opening  dynamics,  but  rather  to  use  test  data  together  with  dimensional  analysis 
techniques  to  develop  correlation  parameters  that  would  inch'de  the  effects  of 
parachute  canopy  stiffness.  However,  to  aid  in  interpreting  the  test  data,  and  to 
guide  in  the  development  of  correlation  parameters,  some  theories  were  examined. 

A  listing  of  some  of  the  early  opening  dynamics  theories  can  be  found  in 
[11-20).  Many  of  these  methods  require  difficult  to  measure  empirical  inputs,  such 
as  the  change  in  projected  canopy  area  with  time.  The  two  theories  that  proved 
most  useful  for  the  current  study  were  those  developed  by  Fu  [5]  and  Lingard  [6]. 
The  dimensional  analysis  work  conducted  by  Fu  was  summarized  earlier  in  this 
report. 


Fu  develops  a  theory  to  predict  parachute  opening  shock  and  time,  using  a 
minimum  of  experimental  inputs  based  on  as  little  test  data  as  possible.  Fu 
considers  the  payload  and  canopy  as  two  point  masses  connected  by  spring  with 
damping.  The  canopy  shape  is  represented  as  a  truncated  cone  topped  Avith  a  hemi¬ 
spherical  cap.  The  apparent  mass  of  the  parachute  is  assumed  to  be  proportional  to 
the  canopy  volume.  The  canopy  vent  size  and  effective  porosity  of  the  canopy  fabric 
are  considered  in  the  conservation  of  mass  equations  used  to  model  canopy  filling. 
Using  Newton's  law,  conservation  of  mass,  and  the  appropriate  geometric  relations, 
a  system  of  non-linear  differential  equations  was  solved  numerically. 

Fu  studied  the  effect  of  the  following  variables  on  opening  time  and  opening 

shock: 

a.  Froude  number 

b.  mass  ratio 

c.  deployment  flight  path  angle 

d.  suspension  line  elasticity  and  damping 

e.  suspension  line  semi-aperture  angle  at  the  beginning  of  inflation 
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All  of  Fu's  comparisons  between  theory  and  experiment  were  for  the  28  ft 
diameter  C-9  parachute.  For  the  test  data  that  Fu  had  available,  satisfactory 
agreement  between  theory  and  experiment  was  obtained.  For  example,  the 
maximum  deviation  of  the  maximum  filling  force  was  less  than  10%  of  the 
measured  values.  This  is  certainly  satisfactory  agreement  in  parachute  technology, 
since  the  variations  occurring  in  experiments  are  usually  considerably  larger. 

Fu's  theory  provides  two  important  graphs  of  relevance  for  the  current  study. 
These  graphs  are  shown  as  Figs.  2  and  3.  Fig.  2  is  a  plot  of  non-dimensional  filling 
time  vs.  mass  ratio  for  various  Froude  numbers  for  vertical  deployments  of  a  C-9 
t\’pe  parachute.  Examination  of  this  figure  shows  that  the  effect  of  Froude  number 
variation  on  non-dimensional  opening  time  is  predicted  to  be  quite  small,  except  at 
very  low  mass  ratios.  Non-dimensional  opening  time  does  vary  with  mass  ratio, 
but  at  mass  ratio  above  0.1,  the  variation  is  small.  As  mass  ratio  increases  from  0.1 
to  10.0,  the  predicted  non-dimensional  opening  time  decreases  only  from  six  to  four. 

As  shown  in  Fig.  3,  the  variation  of  opening  shock  with  both  Froude  number 
and  mass  ratio  can  be  quite  large.  Based  on  these  graphs,  it  may  be  reasonable  to 
look  for  different  correlating  parameters  for  opening  time  than  for  opening  shock. 

Lingard's  analysis  for  predicting  opening  shock  is  considerably  different  from 
that  of  Fu,  but  leads  to  similar  results.  His  method  is  intended  to  provide  a 
relatively  simple  theory  to  let  a  designer  perform  tradeoffs  upon  the  effects  on  the 
peak  load  of  parameters  such  as  snatch  velocity,  suspended  mass,  parachute  size, 
and  altitude  and  angle  of  deployment.  The  main  assumptions  made  by  Lingard  are 
as  follows: 

a.  The  parachute  inflates  in  a  constant  non-dimensional  opening  time 
independent  of  snatch  velocity’  and  mass  ratio. 

b.  For  a  given  parachute  design,  there  exists  a  unique  force  coefficient  vs. 
non-dimensional  time  function,  which  is  independent  of  snatch 
velocity  or  mass  ratio. 

These  two  assumptions  are  then  used  to  extract  empirical  data  from  a  small 
number  of  full-scale  tests  of  a  particular  solid  cloth  parachute  system.  With  these 
assumptions,  Newton's  second  law  is  then  used  to  develop  a  model  that  allows 
calculation  of  force  and  velocity  profiles  during  deployment,  and  the  effects  of  the 
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C-9  parachute  in  vertical  drops. 

(Figure  laKco  (rnm  Rcl.  5) 
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above-mentioned  variables  on  deployment  can  be  studied.  Ungard's  results  for 
non-dimensional  opening  shock  are  shown  in  Fig.  4. 

Comparison  of  Figs.  3  and  4  shows  similar  trends  from  two  totally  different 
approaches  and  lends  confidence  to  the  results  obtained.  Because  of  assumption  (a) 
used  above,  Ungard's  method  is  not  useful  in  predicting  Froude  number  and  mass 
ratio  effects  on  opening  time,  but  the  assumption  itself  agrees  fairly  well  with  the 
theoretical  results  presented  from  Fu  in  Fig.  2. 

In  later  evaluations  of  experimental  data  in  this  report,  the  trends  shown  in 
Figs.  2  through  4  will  be  referred  to  for  guidance  in  evaluating  the  importance  of 
various  parameters  in  data  correlation. 


Effect  of  Canopy  Stiffness  cn  Opening  Dynamics 

The  stiffness  of  the  canopy  affects  the  dynamics  of  the  opening  process  in  at 
least  th’ee  ways:  at  the  time  of  parachute  snatch,  in  the  "lines  taught" 
configuration,  the  canopy  stiffness  affects  the  suspension  line  semi-aperture  angle, 
^0,  and  the  amount  of  open  mouth  diameter,  d,  available  for  the  beginning  of  the 
filling  process  (see  Fig.  1). 

The  stiffness  subsequently  affects  the  "unfolding  process”  of  the  parachute  as 
it  goes  from  the  cylindrical  streaming  configuration  to  a  hemispherical  shape. 
Finally,  the  stiffness  affects  the  amount  of  canopy  and  suspension  line  stretching,  or 
“overexpansion*  that  occurs  once  the  canopy  has  reached  the  hemispherical  shape. 

In  most  parachute  opening  theories,  the  only  way  the  canopy  stiffness  can 
enter  the  theory  is  through  the  initial  value  of  |}o  and  d  at  the  time  of  snatch.  Some 
theories,  for  example  that  of  Payne  (14],  do  consider  suspension  line  stiffness,  but  not 
canopy  stiffness.  In  Ungard's  theory,  the  effect  of  canopy  and  suspension  line 
stiffnesses  are  automatically  included  in  the  theory  because  of  the  way  the 
characteristic  force-time  signature  is  obtained  from  the  test  data.  However,  a 
characteristic  force-time  signature  obtained  for  one  size  canopy-suspension  line 
configuration  may  not  be  valid  for  a  different  size  canopy-suspension  line 
configuration.  Thus,  Ungard's  theory  may  not  work  for  significant  changes  in 
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canopy  or  suspension  line  modulus  of  elasticity,  fabric  areal  density,  or  fabric 
porosity.  In  this  case,  a  theory  more  like  Fu's  or  Payne's  would  be  needed. 


Two  Extremes  of  Stiffness  Effects 

There  are  two  limiting  conditions  in  which  parachute  stiffness  or  flexibility 
would  act  in  different  ways.  One  extreme  would  be  the  case  where  the  parachute  is 
of  large  diameter  and  constructed  from  a  thin  fabric.  In  this  case,  the  parachute 
opening  dynamics  would  be  governed  mostly  by  the  initial  parachute  "open 
mouth"  area  (CdSIi,  plus  the  elasticity  of  the  fabric  itself  as  expressed  by  Young's 
modulus,  E.  In  this  case,  the  concept  of  "flexural  rigidity"  (El  term)  as  used  in  beam 
and  plate  bending  theories  would  have  no  significance. 

The  other  extreme  would  be  the  case  where  the  parachute  is  of  small 
diameter  and  constructed  from  relatively  "thick"  fabric.  In  this  case,  upside  dowrm 
hang  tests  would  show  a  large  initial  parachute  "open  mouth"  diameter,  as 
expressed  in  Heinrich's  original  stiffness  index  formulation.  It  is  expected  that  the 
flexural  rigidity  concepts  used  in  beam  and  plate  theories  might  be  appropriate  here 
in  describing  opening  dynamics.  Test  data  on  the  fabric  crease  resistance  would  be 
useful  here.  This  extreme  is  probably  best  exemplified  in  very  small  wind  tunnel 
model  tests. 

In  each  case  (or  at  each  extreme),  the  question  arises  as  to  which  non- 
dimensional  parameters  can  be  used  to  correlate  opening  dynamics  behavior,  and 
what  can  be  used  as  a  measure  of  the  dividing  line  for  the  two  kinds  of  behavior. 
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DEVELOPMENT  OF  RELATIVE  STIFFNESS  INDEX 


Up  to  this  point  in  time,  it  appears  that  the  only  measure  of  parachute  canopy 
stiffness  that  has  been  used  extensively  is  the  one  developed  by  Heinrich  and 
discussed  earlier  in  this  report  in  the  review  of  Lee's  work  [9].  As  mentioned  earlier, 
the  main  value  of  the  Heiiuich  index  is  to  compare  stiffnesses  of  various  canopies 
and  for  use  as  a  basis  for  developing  model  canopies  that  might  approach  prototype 
canopy  stiffness.  The  Heinrich  index  does  not  allow  one  to  take  data  from  a  stiffer 
model  canopy  and  use  it  to  predict  the  behavior  of  a  more  flexible  prototype,  or  vice 
versa.  In  this  section,  a  new  "relative  stiffness  index"  is  proposed  that  allows  one  to 
perform  these  data  correlations,  especially  for  parachute  opening  time. 

After  considering  available  test  data  and  opening  theories  and  reviewing  the 
dimensional  analyses  conducted  to  date,  U)e  author  concludes  the  main  parameter  that  would 
appear  to  affect  stiffness  differences  between  model  and  full  scale  is  the  "canopy  rigidity" 
or  "relative  stiffness"  term  proposed  by  Kaplun  in  the  form 

(ED 

''■pV2D‘ 


or  in  [8]  as 

_  El  H 
pV2l<  pV2D< 

In  neither  reference  is  any  information  given  on  methods  for  evaluating  canopy 
second  moment  of  area,  I. 

The  approach  used  to  develop  this  concept  in  the  present  work  was  as 
follows: 

The  numerator  is  actually  an  expression  for  the  flexural  rigidity  term  used  in 
beam  bending  theory.  If  a  unit  width  of  parachute  fabric  is  considered  as  a  “heam," 
the  second  moment  of  area  would  be  evaluated  as  shown  in  the  figure  below: 
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t  >  fabric  thickness 


I 


1 


—  1 


1 

T 


(l)t3 

12 

Beam  theory,  of  course,  does  not  apply  unless  t  >  width.  Since  the  second 
moment  of  area  I  (commonly  called  moment  of  inertia),  dimensionally  must  have 
units  of  length  to  the  fourth  power,  the  expression  for  fl  was  modified  by  writing  it 
in  the  form 


pV2D< 


Et3 

12  pV2d3 


The  characteristic  velocity  was  selected  as  the  snatch  velocity,  V$,  and  the  numerical 
factor  of  12  was  dropped  (a  valid  step  in  dimensional  analysis). 


The  expression  then  becomes 


(1) 


Furthermore,  it  was  felt  that  canopy  stiffness  would  be  more  analogous  to  plate 
bending  stiffness  than  to  beam  bending  stiffness,  since  the  Poisson's  ratio  effect  is 
neglected  in  the  bending  of  beams,  but  is  important  in  the  bending  of  plates  where 
lateral  contraction  is  inhibited  much  more.  Hence,  the  parameter  E  (Young's 
modulus)  was  replaced  by  the  term  E/(l*v2)  giving  the  final  form  for  the  expression 


E 

pv|  (l-u2) 


(2) 


This  parameter  is  then  similar  dimensionally  to  the  expression  often  used  in 
plate  theory,  i.e.,  the  work  of  Calligeros  (21].  Certainly,  a  parachute  canopy  does  not 
exhibit  the  stiffness  exhibited  by  a  plate  in  small  deflection  theory;  nevertheless. 
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there  must  be  some  similarity  between  the  unfolding  of  fabric  and  the  bending  of  a 
plate.  Furthermore,  analysis  shows  that  if  only  the  overexpansion  of  the  inflated 
canopy  were  considered,  then  the  theory  for  strain  in  a  hemispherical  shell  could  be 
applied,  and  an  expression  of  the  form 


c*  = 


E 

pVj  (1-u) 


(3) 


would  result.  Thus  the  difference  between  Eqs.  (1),  (2)  and  (3)  is  in  the  relative  effect 
of  Poisson's  ratio  on  canopy  rigidity. 

The  form  of  the  parameter  selected  for  study  was 


This  parameter  will  be  called  the  “relative  stiffness  index." 


Modification  of  Form  for  Relative  Stiffness  Index 

The  relative  stiffness  index  given  in  Eq.  (2)  contains  canopy  design  variables 
(E,u,  t.  Dp)  and  drop  condition  variables  (p,Vj).  More  insight  into  the  effect  of  each 
parameter  may  be  obtained  by  modifying  the  expression  as  follows:  If  the  Froude 
number  is  introduced,  the  expression  can  be  changed  from 

''  pv^  W 


using  Fr 


gDoFr  into 


r 


E 

(l-v2) 


f  ’  1 

't3  ' 

iPgPrj 

^  / 
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Multiplying  the  numerator  and  denominator  by  the  parachute  canopy  weight,  Wp, 
and  rearranging  gives 


r__J _ 

^  Wn(l-u2)  iFrJ  W 


Wp  (l-u2) 


gpDo 


mr 


Since  Wp/g  =  mp,  the  parachute  canopy  mass,  and  defining  M^P  =  — 3  as 

pDo 


done  by  Lee  [9],  gives 


(4) 


Wp  2 

Using  Wc  =  — 2 —  average  areal  density  of  the  canopy  gives  Wp  =  JiDoWc/4. 


(nD;/4) 


Substituting  this  into  Eq.  (4)  gives 


rjLVMLE 

Fr 


\  j 

Dropping  the  A/n  term  since  it  is  just  a  constant  like  the  value  of  12  from  plate 
bending  theory  gives  the  final  form 


The  term  in  brackets  is  a  function  only  of  the  parachute  construction,  while 
the  expression  in  parentheses  includes  the  effect  of  drop  conditions,  as  well  as 
canopy  mass  and  size. 
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We  can  define  a  modified  canopy  stiffness  index,  call  it  T]',  different  from 
Heinrich's  stiffness  index+  (but  related  to  it)  as  follows: 


^  "wc  [(i.v2)JiDoJ 

With  this  definition,  the  relative  stiffness  index  parameter  becomes 


(5) 


(6) 


Physical  Significance  of  Relative  Stiffness  Index 

Examination  of  Eq.  (6)  for  the  relative  stiffness  index  term  will  show  the 
importance  of  each  variable  physically.  This  equation  shows  how  the  relative 
stiffness  index  fits  in  with  the  complete  picture  of  relative  canopy  rigidity.  The 
Froude  number  term  in  the  denomination  makes  the  "absolute"  or  "geometric" 
stiffness  into  a  "relative"  stiffness,  since  it  is  now  relative  to  inertia  and  gravity 
forces  during  opening.  The  term  in  the  numerator  Mr,p  is  proportional  to  the  ratio 
of  canopy  mass  to  included  mass:  mp/pD^.  The  importance  of  this  term  was  also 
recognized  by  Kaplun  [1]  as  early  as  1951.  He  called  this  term  the  "canopy  inertia 
number,"  I,  defined  as 


1 


m 

pD3 


where  m  was  the  mass  of  the  parachute  canopy  (Kaplun's  study  dealt  with  "Infinite 
mass"  inflations  only,  thus  the  payload  mass  did  not  enter  the  study  as  a  variable). 
Kaplun  felt  that  it  was 


’  Heinrich’s  stiffness  index  w*s  defined  as »)  ^  Since  Dinaxcf(E,v,t,De),  the  stiffness 

indices  are  obviously  related  in  some  physical  way,  but  in  a  way  probably  not  determinable  by  any 
simple  mathenrutical  analysis. 
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...particularly  important  to  establish  the  effect  of  the  inertia  number,  since  an  opening 
governed  by  canopy  inertia  forces  would  be  characteristically  different  from  one  in  which  the 
mass  of  the  canopy  has  no  effect.  In  the  former  case,  the  rate  of  opening  would  be  determined  by 
the  dynamic  pressure  and  the  inertia  of  the  canop>y,  while  in  the  latter  case  by  the  volume  in 
flow  of  air.  It  is  the  latter,  flow-determined  case  which  has  been  assumed  by  O'Hara  and  by 
W.  Muller  in  the  rational  analyses  which  they  respectively  suggested. 

Thus,  the  relative  stiffness  index  of  Eq.  (2)  includes  the  effect  of  canopy  inertia 
as  well.  To  the  author's  knowledge,  the  present  study  is  the  first  one  to  identify  the 
importance  of  the  variables  in  the  combination  shown  in  Eq.  (6). 


Overview 

The  next  step  of  the  present  study,  then,  is  to  examine  available  parachute  test 
data  to  see  if  stiffness  effects  can  be  correlated  using  the  relative  stiffness  index  of  Eq. 
(2).  The  study  will  be  based  primarily  on  the  exp)erimental  data  of  Lee  (9,22),  as  this  is 
the  only  data  to  the  author's  knowledge  where  several  different  fabric  densities  were 
tested  on  several  parachute  scales  as  well.  A  small  amount  of  wind  tunnel  test  data 
will  also  be  examined,  as  well  as  some  full-scale  drop  tests,  which  were  conducted 
using  similar  parachutes  constructed  with  fabrics  of  different  areal  densities. 

In  order  to  conduct  the  desired  data  correlation  studies,  it  is  first  necessary  to 
develop  data  for  the  parachute  fabric  properties  needed  in  the  relative  stiffness 
index.  This  task  involves  obtaining  data  for  Young's  modulus  and  Poisson's  ratio  for 
various  parachute  fabrics  and  will  be  examined  in  the  next  section  of  this  report. 
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EVALUATION  OF  FABRIC  PROPERTIES 


To  use  the  relative  stiffness  index  it  is  necessary  to  have  values  of  Young's 
modulus,  E,  and  Poisson's  ratio,  u,  for  parachute  fabric.  To  truly  characterize  the 
properties  of  parachute  fabric,  it  must  be  considered  as  a  nonlinear  material  with 
complex  strain,  strain  rate,  apd  hysteresis  properties.  If  the  fabric  is  assumed  to 
behave  linearly,  the  general  stress-strain  relationships  for  any  linearly  elastic, 
anisotropic  material  must  be  represented  by  six  equations  which  relate  the  six  most 
general  components  of  strain  to  the  six  components  of  stress. 

Somewhat  less  complicated  representations  of  fabric  properties  have  been 
conducted  by  Steeves  [23],  and  by  Remington,  Madden,  and  O'Callahan  [24]. 
However,  these  models  were  for  specialized  applications,  and  still  too  complex  to 
use  for  this  study  of  scaling  parachute  stiffness.  Since  fabrics  are  not  isotropic 
materials,  strictly  speaking,  there  is  no  single  value  of  E  and  u  that  can  be  used  to 
characterize  fabric  properties.  However,  what  is  needed  is  some  tj'pe  of  estimate  for 
parachute  cloth  in  order  to  apply  the  method. 

Fabric  data  from  several  sources  were  examined  and  estimates  for  E  and  u 
attempted.  Since  parachute  fabric  is  stressed  biaxially  during  canopy  opening  and 
descent;  it  would  be  desirable  ideally  to  obtain  these  values  from  biaxial  test  data  for 
each  particular  fabric  of  interest.  Both  biaxial  and  uniaxial  test  data  were  examined 
and  isotropic  and  orthotropic  models  for  fabric  properties  were  each  evaluated. 


Biaxial  Test  Data 

The  most  consistent  results  for  both  Young's  modulus,  E,  and  Poisson's  ratio, 
u,  were  obtained  from  biaxial  test  data  obtained  for  50  g/m^  (1.47  oz/yd-)  plain 
weave  British  nylon  parachute  fabric,  and  reported  by  Render  and  Bradley  [25]. 

The  fabric  was  assumed  to  behave  isotropically,  in  which  case  the  strain 
equations  in  the  two  directions  would  be  given  by 
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ei  =£(oi  -“002) 

C2=f  (02*‘D0l) 

These  equations  were  solved  for  E  and  v  to'give 

2  2 
Oj  -02 

E  = - 

OlEl  -  02E2 

02E1  *  OlE2 

\)  = - 

OlEl - 02E2 

Eight  test  points  were  chosen  at  random  from  the  biaxial  stress-strain  data  as 
reported  by  Render  and  Bradley  and  these  points  are  shown  on  their  plot  in  Fig.  5. 
Notice  that  some  data  points  at  negative  values  of  strain  were  also  taken. 

The  data  were  substituted  into  Eqs.  (7)  and  (8)  and  the  results  obtained  are 
shown  in  Table  1.  It  can  be  seen  that  the  data  are  fairly  consistent,  and  the  standard 
deviations  obtained  are  reasonable  under  the  circumstances.  This  led  to  confidence 
that  this  method  would  work  for  other  fabrics. 

Table  1.  Properties  of  1.47  oz/yd^  plain  weave,  British  nylon  parachute  fabric 


(7) 

(8) 


Point  Selected 
(sec  Fig  5) 

Young's  Modulus 

E  N/mm 

Poisson's  Ratio 

V 

1 

38.5 

0.154 

2 

33.9 

0.228 

3 

35.5 

0.346 

4 

40.5 

0.215 

5 

41.7 

0.095 

6 

39.3 

0.130 

7 

41.5 

0.225 

8 

36.4 

0.273 

Arrange  values: 

38.4  N/nvn 

0.208 

Standard  deviation: 

2.9  N/mm 

0.081 
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A  similar  approach  was  tried  on  the  biaxial  stress  strain  data  for  four  fabrics 
tested  by  Monego  and  Pilsworth  [26]  using  a  pressurized  cloth  cylinder  for  the  tests. 
It  was  hoped  that  this  technique  would  yield  values  for  heavier  fabrics  that  could  be 
used  to  estimate  Poisson's  ratio  for  the  heavy  fabrics  used  for  the  parachute  drop 
tests  reported  by  Lee.  Unfortunately,  the  isotropic  stress  strain  equations  would  not 
give  meaningful  values  of  Young's  modulus  or  Poisson's  ratio  for  any  of  the 
heavier  fabrics.  It  was  next  decided  to  look  at  uniaxial  test  data. 


Uniaxial  Fabric  Test  Data 

From  the  [26]  test  data  for  uniaxial  tests  to  failure  on  conventional  weight 
parachute  fabric,  it  was  possible  to  obtain  separate  values  for  stiffness  in  the  warp 
and  fill  directions.  Assuming  a  fabric  thickness  of  0.0028  in  as  measured  in  Lee's 
later  tests,  the  following  secant  modulus  values  were  obtained: 

Ewtrp  “  7.80  X  10^  psi 

Efiii  *  5.81  X  ICH  psi 

Have  *  6.8  X  1  O'*  psi 

An  order  of  magnitude  for  Poisson's  ratio  for  this  fabric  was  assumed  by  modifying 
the  isotropic  equations  to  the  following  form: 


Ow  ^ 

^''■’Ew’  Ef 


(9) 


VOw  Of 

Ew  ■  Ef 


(10) 


The  use  of  two  different  values  of  Young's  modulus  in  this  form  of  equation 
violates  the  isotropic  assumptions  on  which  the  equations  are  based.  However,  use 


O  w  1  O  2 

of  these  equations  on  four  test  data  points  at  stress  ratios  of  —  ■  r  and  —  «  7  gave 

Of  *  Of  * 

values  of  Poisson's  ratio  in  the  range  from  0.045  <  u  S  0.270,  with  an  average  value 
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of  u  *=  0.14.  This  value  compares  favorably  with  the  value  of  0.21  for  the  British 
parachute  fabric. 


Linear,  Orthotropic  Fabric  Model 

Because  the  above  technique  yielded  some  promising  values  for  Poisson's 
ratio,  it  was  decided  to  model  the  fabric  more  accurately  using  the  following  linear 
orthotropic  stress-strain  relations: 

Ex  =  SnOx +  Si20y  (11) 

£y  =  S2iOx  +  S220y  02) 

In  this  model,  Sn  and  S22  can  be  different  from  each  other,  but  it  is  required  that 
S12  '  S21.  The  isotropic  equations,  arranged  in  similar  form,  are 


Comparing  these  equations  with  Eqs.  (11)  and  (12)  shows  that  the  orthotropic  model 
allows  the  use  of  separate  values  of  Young's  modulus  in  the  warp  and  fill  directions, 

but  the  requirement  that  S12  =  S21  would  require  the  use  of  a  term  like  for 

$12  and  S21. 


The  values  of  Sn  and  S22  would  come  from  uniaxial  fabric  tests  in  the  warp 
and  fill  directions  and  we  would  then  try  to  represent  the  fabric  properties  as 


Using  the  values  of  Ew,  Ef,  and  E*ve  from  the  uniaxial  test  data,  the  equations 
given  above  were  applied  to  both  the  pressurized  cylinder  and  cruciform  test  data 
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for  the  MIL-C-7020-1  parachute  fabric  and  the  results  are  given  in  Table  2.  Note  that 
the  value  of  0.142  from  Table  2  compares  favorably  with  the  value  of  0.14  reported 
earlier  using  Eqs.  (9)  and  (10). 

Table  1  Values  of  Poisson's  ratio  for  MIL-C-7020-1  type  parachute  fabric 
using  Eqs.  (131  and  (14) 


Test  Type 

Poitton't  Ratio,  v 

0%^ 

Cruiform,  —  «  1.0,  Warp 

Of 

0f)880 

Ow 

Cniiform,  —  *1.0,  Fill 

Of 

0.234 

Press.  Cvl,—  ■  1.0,  Warp 

•  Of 

0.170 

Press.  Cvl,  —  .  1.0,  Fill 

'  Of 

0.0638 

Cruiform,  —  ■  2.0,  Warp 

0.210 

Cruiform,  —  »  2.0,  Fill 

Of 

0.038 

Press.  Cvl,  —  ■  2.0,  Warp 

0.207 

Press.  Cyl,  —  ■  2.0,  Fill 

0.129 

Avcracc  Value: 

0.142 

An  attempt  to  use  this  approach  (modified  orthotropic  equations)  for  the  body 
armor  fabric  of  (26)  was  tried,  but  the  values  obtained  for  Poisson's  ratio  seemed 
unreasonable. 

The  final  choice  for  values  of  Young's  modulus  was  based  on  uniaxial  test 
data  obtained  for  samples  of  the  cloth  used  to  construct  the  parachutes  used  in  each 
of  Lee's  drop  tests.  Average  values  from  tests  in  the  warp  and  fill  directions  were 
used  for  each  of  the  fabric  weights  given  in  Table  3  below. 
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Table  3.  Uniaxial  test  data  for  Young's  Modulus  and  average  values  used  for  analysis 


The  final  input  values  to  the  relative  stiffness  index  were  to  use  the  values  of 
average  Young's  modulus  from  the  table  above  for  each  individual  fabric,  and  a 
common  value  of  Poisson's  ratio  u  =  0.14  for  all  fabrics. 

Ordinary  changes  to  the  value  of  Poisson's  ratio  for  the  different  fabrics,  if  the 
values  could  be  found,  would  have  negligible  effect  on  the  final  numerical  value  of 
relative  stiffness  index.  The  scatter  in  experimental  data  would  obscure  any  attempt 
to  verify  the  Poisson's  ratio  effect  accurately. 
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APPLICATION  OF  MEITiOD  TO  EXPERIMENTAL  DATA  CORRELATION 


C-9  Parachute  Test  Data 

The  most  systematic  experimental  study  of  the  effects  of  fabric  stiffness  on 
parachutes  of  various  sizes  is  the  work  done  by  Lee  (9,22).  In  [9],  Lee  reports  the 
results  of  drop  tests  conducted  on  full-scale  28  ft  diameter,  solid  cloth  C-9  parachutes, 
as  well  as  for  a  14  ft  diameter  l /2-scale  model  and  a  7  ft  diameter  V4-scale  model. 
The  three  parachutes  were  constructed  of  1.1  oz/yd^  ripstop  nylon. 

In  [22],  Lee  presents  similar  drop  test  data  on  full-scale  and  V4-scale  C-9 
parachutes  using  fabric  of  the  following  weights:  1.1  oz,  4.1  oz,  9.2  oz,  and  12.2 
oz/yd2.  Lee  presents  data  for  non-dimensional  opening  force  and  non-dimensional 
opening  time  vs.  Froude  number,  Fr,  for  two  mass  ratios. 

A  matrix  of  the  test  conditions  is  presented  in  Table  4. 

Table  4.  Range  of  test  vatiables  covered  in  Lee's  tests 


Scale 

Actual  Fabric  Areal 
Densttv,  oz/vd^ 

Froude  Number 
Range 

M,,i  ■  0.U  (Ref.  22) 

Full  (C-9) 

0.94 

4—7 

Full  (C-9) 

4.1 

4.5—5 

Full  (C-9) 

9.2 

2—7.5 

Full  (C-9) 

12.2 

2.5-4 .5 

One-Quarter 

0.94 

3.5-50 

One-Quarter 

4.1 

4—20 

Mrj  -  0.20  (Ref.  9) 

Full  (C-9) 

0.94 

4—6 

One-Half 

0.94 

4—20 

One-Quarter 

0.94 

4—60 

All  drop  tests  were  done  outside,  from  a  crane  for  vertical  drops,  and  from  a 
packed  configuration.  Lee's  test  data  are  presented  for  operung  time,  as  opposed  to 
filling  time.  Opening  time  is  here  defined  as  the  time  from  occurrence  of  first 
snatch  force  to  the  time  of  peak  opening  force  (filling  time  is  usually  defined  as  the 
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time  at  which  the  canopy  first  reaches  a  projected  area  equal  to  the  steady  state 
projected  area). 


Opening  Time  Data  Correlation,  Mass  Ratio  0.13 

Fig.  6  is  a  composite  plot  of  all  the  data  presented  in  [22]  for  non-dimensional 
opening  time,  t,  vs.  Froude  number.  It  can  be  seen  that  the  correlation  shown  is 
fair.  The  effects  of  fabric  areal  density  relative  to  canopy  diameter  are  somewhat 
difficult  to  sort  out,  but  some  conclusions  can  be  reached. 

The  data  points  for  the  smallest  parachute  with  a  relatively  heavy  fabric  0/4 
scale,  4.1  oz/yd^)  have  the  shortest  non-dimensional  opening  time,  falling  below 
nearly  all  the  other  data.  The  data  for  the  full-scale  parachute  using  the  lightest  (1.1 
oz/yd^)  fabric  (and  even  with  the  4.1  oz/yd^  fabric)  show  the  longest  opening  time. 
The  V4-scale  parachute  with  the  light  fabric  and  the  full-scale  parachutes  with  the 
heavy  fabrics  fall  somewhere  in  the  middle. 

If  all  these  data  points  are  plotted  against  the  relative  stiffness  index,  the 
correlation  shown  in  Fig.  7  results.  It  can  be  seen  that  this  correlation  is  considerably 
better  than  that  shown  in  Fig.  6. 


Opening  Time  Data  Correlation,  Mass  Ratio  0.20 

Examination  of  Fig.  2  shows  that  a  change  in  mass  ratio  from  0.13  to  0.20 
would  have  only  a  very  small  effect  on  theoretical  filling  time,  much  less  than  the 
range  of  experimental  data  scatter  for  a  given  mass  ratio.  It  is  expected  that  the  effect 
on  opening  time  would  be  just  as  small. 

With  ,;us  logic  in  mind,  the  data  from  [9]  for  a  mass  ratio  of  0.20  were  plotted 
together  with  the  mass  ratio  0.13  data.  This  provides  additional  data  points  for  the 
V4-scale  and  full-scale  C-9  parachutes,  as  well  as  data  for  the  V2-scale  C-9  parachute. 
The  correlation  for  non-dimensional  opening  time  vs.  relative  stiffness  index  is 
shown  in  Fig.  8  (the  flagged  data  points  are  for  mass  ratio  0.20).  The  correlation  is 
still  seen  to  be  generally  good.  By  contrast,  if  the  other  mass  ratio  is  simply  plotted 
vs  Froude  number,  the  graph  in  Fig.  9  results. 
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1  2  34Sf7t  810  20  to 

Vs* 

Fro\jde  number,  Fr  ■ 

gDo 

Fig.  6;  Opcrung  Dmc  as  a  funaion  of  Froude  number  for  full-scale  and  ^  /4-scale  C-9  parachutes  at  mass  rabo  0.13. 
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function  of  relalivc  stiffness  index  for  full-scale  and  >/l-5caleC-9paracl»ulcs  at  mass  ratio  013  Rt-I.itive  stiffness 


Kcl.ilivc  stiffness  iiulex. 


9  4  567a910 


=t'-T- 


20  90  40  SO  00 


Vs* 

Fjg.  9;  Opening  time  as  a  funcoon  of  Froude  number  (or  full-scale,  proude  nuirto,  Fr  ■ 

^  /2-scale,  ar»d  ’  /«-$calc  C-9  parachutes  ai  mass  ratios  0.13  and  0.20.  ®  ® 
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In  Fig.  8,  there  are  a  few  data  points  that  correlate  poorly  in  the  range  C  =  I  I 
1.8  X  ID’S.  The  reason  for  this  is  not  known.  Note,  however,  that  the  V2*scale  data 
fit  nicely  between  the  V4-scale  and  full-scale  data.  The  main  interesting  trend 
shown  here  experimentally  agrees  with  Fu's  theory  in  that  there  is  no  discernable 
difference  between  the  two  mass  ratios,  and  that  the  effects  of  parachute  scale  and 
fabric  stiffness  are  greater  than  the  effects  of  change  in  mass  ratio. 


Wind  Tunnel  Data 

In  [20],  Heinrich  presents  test  data  for  small-scale  wind  tunnel  models  of  the 
C-9  parachute  that  were  3  ft  in  diameter.  A  finite  mass  test  was  simulated  by  using  a 
suspended  mass  attached  to  a  cable.  As  the  parachute  opened  in  the  wind  tunnel,  it 
moved  downstream  in  the  tunnel,  simultaneously  lifting  the  mass.  The  mass  ratio 
used  was  0.24,  again  close  enough  to  the  values  of  0.13  and  0.20  previously  plotted, 
considering  experimental  data  scatter. 

Heinrich  measured  filling  time  as  the  point  where  projected  canopy  area  first 
equals  steady  state  projected  area.  The  data  points  are  shown  on  Fig.  8.  The  points 
fall  just  above  the  data  points  for  the  4.1  oz/yd^,  V4-  scale  C-9  canopy,  but  if  the 
opening  time  were  known  as  opposed  to  the  filling  time,  the  points  would  be  even 
closer  to  the  other  data  points.  The  difference  between  these  values  and  the  other 
values  may  also  be  due  to  two  other  reasons:  the  manner  in  which  the  wind  tunnel 
test  is  conducted  makes  it  difficult  to  simulate  a  vertical  drop  test,  and  at  this  small  a 
scale  factor  (approximately  VlO),  the  effect  of  reinforcing  along  the  parachute  top 
cord  lines  may  overshadow  the  basic  stiffness  characteristics  of  larger,  solid  flat 
circular  parachutes. 

It  is  interesting  to  note  that  all  the  data  points  on  Fig.  8  are  correlated  fairly 
well  by  the  equation; 

T  «  0.86(;)-O  064  (15) 

This  equation  covers  the  complete  range  of  mass  ratios  and  Froude  numbers 
considered  in  all  C-9  tests  discussed  to  this  point. 
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Application  of  the  Method  to  35  ft  Do  T-10  Parachutes 

Watson  [27]  presents  test  data  for  a  series  of  drop  tests  for  the  35  ft  diameter  T- 
10  parachute.  A  series  of  46  or  more  drops  were  conducted  from  airaaft  at  altitudes 
ranging  from  approximately  800  ft  to  38,000  ft;  snatch  velocities  varying  between  200 
and  475  fps;  and  suspended  loads  ranging  from  174  to  865  lb.  Two  typ>es  of  canopy 
construction  were  used:  a  standard  1.1  oz/yd^  (nominal)  T-10  canopy,  and  a  high- 
strength  (2.25  oz/yd2)  extended  skirt  canopy.  Various  trajectory  deployment  angles 
were  also  tested. 

Watson's  data  were  not  presented  in  dimensionless  form,  so  it  was  necessary 
to  convert  the  data  into  the  format  desired.  In  addition,  the  time  data  were  not 
presented  as  an  opening  time  (time  from  snatch  force  to  peak  force),  but  as  a  filling 
time  (time  to  first  full-open  projected  area).  Therefore,  a  non-dimensional  filling 
time,  “tf,  was  plotted  vs.  the  stiffness  parameter, 

The  data  are  presented  in  Fig.  10,  together  with  one  data  point  obtained  on  a 
V  10-scale  model  wind  tunnel  test  conducted  by  Heinrich  [28]  at  finite  mass 
conditions.  The  data  correlation  in  this  case  does  not  follow  the  same  trend  as 
shown  in  Fig.  8,  that  is,  a  decreasing  value  of  t  with  increasing  relative  stiffness 
index,  (j.  As  a  matter  of  fact,  except  for  the  single  point  for  the  wind  tunnel  test,  it 
appears  that  the  non-dimensional  filling  time  is  essentially  constant — independent 
of  stiffness  index,  payload  mass  ratio,  or  trajectory  inclination  angle. 

To  explore  this  phenomenon  further,  Watson's  original  data  were  replotted 

several  different  ways.  Watson's  original  plot  of  filling  time  vs.  snatch  velocity  is 

reproduced  as  Fig.  11.  Here  it  can  be  seen  that  filling  time  decreases  as  snatch 

velocity  increases.  It  can  be  seen  that  the  high  strength  canopy  has  a  slightly  longer 

filling  time  than  the  conventional  T-10  canopy.  When  the  same  data  are  non- 

dimensionalized,  a  plot  of  non-dimensional  filling  time,  tf,  vs.  Froude  number  Fr  = 
2 

Vj/gDo  is  obtained,  as  shown  in  Fig.  12.  It  now  appears  that  the  parachute  inflations 
show  the  constant  (non-dimensional)  filling  time  characteristic  pointed  out  by 
earlier  investigators. 

To  further  sort  out  any  trends,  the  data  were  replotted  once  more  in  Fig.  13, 
where  non-dimensional  filling  time  is  plotted  vs.  "Froude  number  of  operation," 


35 


36 


filling  liiiK*  as  a  fwiKlion  of  relative  stiffness  index  for  TdO  and  other  parachute  configurations 


Filling  Time  (tf).  Seconds 


200  400  600 


Snatch  Velocity  IVj),  ft/sec 

rig.  11:  Filling  timf  vs.  sn»tch  velocity  for  T-IO  extended  skin  canopv 
and  high  strength  canopy. 

35  ft.  Do  extended  skirt  canopy  and  high  strength  canopy. 
(Figure  taken  from  Ref.  27) 


37 


T-10  CANOPY 

O  174  lb  &  193  lb 
^  207  lb 

250  lb  Horiz  Traj 
(t>  250  lb  Ven  Traj. 


H.S,  CANOPY 
♦-250  lb  HoriiTtaj. 
^  250  lb  Vert  Traj. 
^  480  lb  Hori:  Traj. 
•■865  lb  Horiz  Traj, 


■*^480  lb  Horiz  Traj.  $  865  lb  Vert  Traj.' 


3  4  3  i  7  •  {  I  2  3  4 

F, 

Fig.  12;  Non-dimensional  filling  time  v$.  Froude  number  for 
T-10  extended  skin  cartopy  and  high  strength  canopy. 
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Fr/siny.  In  this  manner,  the  horizontal  scale  is  expanded,  but  no  new  trends  are 
observed.  The  flight  path  angles  at  deployment  ranged  from  "horizontal 
deployments"  with  trajectory  snatch  angles  as  flat  as  4®  below  horizontal,  to  "vertical 
deployments"  with  trajectory  snatch  angles  as  steep  as  75®  below  horizontal.  The 
filling  time  appears  to  be  independent  of  snatch  angle  when  non-dimensionalized 
vs.  Froude  number  of  operation. 

Payload  mass  ratios  for  these  tests  varied  from  as  low  as  0.063  to  as  high  as 
0.70.  The  non-dimensional  opening  time  being  independent  of  mass  ratio  agrees 
with  Fu's  theoretical  results  for  the  C-9  parachute  shown  in  Fig.  2.  Here  filling  time 
is  essentially  independent  of  mass  ratio  for  mass  ratios  above  approximately  0.1.  It  is 
expected  that  the  T-10  parachute  would  show  the  same  general  characteristics. 

The  only  parameter  that  seemed  to  affect  filling  time  in  a  statistically 
significant  way  was  the  fabric  density.  The  shaded  symbols  for  the  heavier  (2.25 
oz/yd2)  construction  show  a  longer  filling  time  than  the  standard  1.1  oz/yd2  fabric. 
The  average  non-dimensional  filling  times  are  y  *  8.4  for  the  standard  construction 
parachute,  based  on  28  tests;  and  tf  »  9.6  for  the  heavier  construction  parachute, 
based  on  18  tests.  These  values  are  shown  on  Fig.  13.  Note,  however,  that  the  data 
scatter  within  either  set  of  fabric  tests  is  greater  than  the  difference  between  the  two 
mean  values. 

Going  back  to  the  T-10  data  correlation  shown  in  Fig.  10  and  comparing  it 
with  the  C-9  correlation  shown  in  Fig.  8,  it  appears  that  there  was  not  enough 
variation  in  the  relative  stiffness  index  to  show  any  trend  for  the  T-10  parachute. 
This  was  primarily  because  the  variation  in  parachute  and  test  conditions  did  not 
cover  a  wide  enough  range.  The  thickness  to  diameter  ratio  of  the  parachutes 
varied  only  from  6.7  x  lO*^  to  1.6  x  lO'^,  as  opjX)sed  to  from  8.3  x  10"^  to  1.3  x  10'^  for 
Lee's  C-9  ty’pe  drop  tests. 

It  is  interesting  to  note  that  passing  a  line  through  the  T-10  wind  tunnel  data 
point  and  using  the  same  slope  as  the  line  on  the  C-9  data  plot  correlates  the  T-10 
tunnel  test  point  well  with  the  full-scale  data.  This  line  is  shown  in  Fig.  10. 
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Large  Diameter  Parachutes  (G*12  and  G*ll) 

Test  data  were  available  to  plot  single  points  for  the  64  ft  D©  G-12D  parachute 
and  the  100  ft  Do  G-11  parachute.  The  available  data  were  incomplete  for  each  of 
these  parachute  drops,  and  specification  values  had  to  be  used  for  fabric  properties, 
rather  than  actual  values  from  measurements.  Filling  times  rather  than  opening 
times  were  non-dimensionalized  and  plotted  in  Fig.  10.  Drop  altitudes  also  had  to 
be  estimated. 

The  G-12D  data  point  (2.25  oz/yd^  fabric)  falls  midway  between  the  two  groups 
of  data  shown  for  the  T-10  parachutes.  The  payload  mass  ratio  for  th  s  drop  was  0.13, 
within  the  range  of  mass  ratios  covered  by  the  T-10  drop  tests.  For  the  G-12D 
diameter  and  fabric  weight  used,  the  performance  of  this  parachute  puts  it  into  the 
area  of  constant  non-dimensional  filling  time,  based  on  one  data  point. 

The  data  point  for  filling  time  of  the  100  ft  D©  G-11  parachute  was  for  a 
payload  mass  ratio  of  0.06.  This  is  somewhat  lower  than  the  lower  mass  ratios  tested 
for  the  T-10  parachute,  but  the  data  point  is  seen  to  fall  along  an  extension  of  the 
line  connecting  filling  time  for  the  T-10  wind  tunnel  model  with  filling  time  for  the 
heavier  fabric  T-10  prototypes. 
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OPENING  SHOCK  CORRELATION  PARAMETERS 


The  approach  itviiialiy  taken  for  analyzing  parachute  opening  shock  during 
filling  was  the  same  as  that  used  to  analyze  opening  time.  Available  theories  were 
used  as  a  guide  to  interpreting  test  data.  The  scatter  in  test  data  for  nominally 
identical  drop  conditions  again  makes  it  difficult  to  sort  out  trends  unless  test 
parameters  are  changed  significantly,  almost  drastically. 

A  review  of  Figs.  3  and  4  from  theory  show  that  mass  ratio  and  Froude 
number  have  a  significant  effect  on  non-dimensional  opening  shock.  Fig.  2  shows 
that  these  parameters  have  a  small  effect  on  opening  time.  For  this  reason,  it  ntight 
be  expected  that  the  correlation  parameters  for  stiffness  effects  might  affect  opening 
shock  differently  than  for  opening  time. 


Correlation  Study  for  Opening  Shock 

Several  different  approaches  were  tried  to  correlate  opening  shock  force  to 
parachute  canopy  stiffness  Most  of  these  attempts  were  unsuccessful.  The  first 
attempt  was  to  correlate  non-dimensional  opening  shock  directly  with  the  relative 
stiffness  index,  Since  opening  time  correlated  well  with  relative  stiffness  index,  it 
was  expected  that  opening  shock  would  correlate  well  also.  It  turned  out  that  the 
correlation  was  unsatisfactory. 

A  number  of  other  attempts  were  tried  for  various  correlations  using 
modifications  of  the  relative  stiffness  index,  for  example,  raising  the  ratio  of  fabric 
thickness  to  parachute  diameter,  t/D©,  to  the  first  or  second  power,  instead  of  to  the 
third  power.  None  of  these  attempts  were  successful  either. 

Since  the  peak  opening  load  occurs  during  the  last  portion  of  the  opening  time  interval, 
it  may  be  affected  differently  by  parachute  line  elasticity  and  canopy  elasticity  than  opening 
time.  Much  of  the  time  during  inflation  the  forces  in  the  suspension  lines  and  canopy  are 
relatively  km*.  It  is  only  during  the  final  monroents  of  inflation  that  the  suspension  line  and 
canopy  forces  become  large.  Ths  fact  would  explain  the  difference  in  the  way  mass  ratio  and 
Froude  number  affect  opening  shock  as  opposed  to  affecting  opening  time. 


42 


Effect  of  Sxispension  Line  Elasticity 

Goodrick  [29]  made  a  simplified  study  of  the  effects  of  installing  a  spring  and 
damper  between  a  parachute  and  its  payload.  This  study  assumed  that  the  parachute 
was  already  initially  fully  inflated  at  the  snatch  velocity.  He  studied  the  effect  of 
numerical  variation  in  the  spring  constant  and  damper  coefficient  on  peak  opening 
shock  load.  In  this  man:  3r,  the  filling  time,  Xf,  was  set  at  zero,  but  different  values 
of  opening  shock  occurred  for  varying  payload  mass,  snatch  velocity  and  simulated 
suspension  characteristics.  It  was  found  that  the  shock  experienced  by  the  payload 
was  affected  greatly  by  the  simulated  suspension  characteristics. 

In  Payne's  theoretical  study  [14],  suspension  line  elasticity  was  varied  together 
with  initial  suspension  line  stretch,  and  it  was  found  that  the  canopy  opening 
characteristics  were  affected  by  these  values.  Payne's  theory  was  the  first  one  to 
predict  the  characteristic  two  peaks  in  the  force-time  history.  It  was  found  that  riser 
and  shroud  line  stiffness  and  initial  stretch  have  a  significant  effect  on  opening 
history  and  loads. 


Modified  Impulse  Concept 

The  basic  physics  of  the  canopy  opening  process  and  payload  deceleration 
dictate  the  following.  The  payload  mass,  snatch  velocity,  and  terminal  velocity  are 
the  only  quantities  that  determine  the  total  retarding  impulse  that  must  be  applied 
to  the  payload  to  decelerate  it  to  the  desired  velocity.  Neglecting  payload  drag,  the 
required  impulse  can  be  achieved  by  a  large  retarding  force  occurring  over  a  short 
period  of  time,  or  a  smaller  retarding  force  occurring  over  a  longer  period  of  time. 
This  principle,  applied  to  the  payload  itself,  is  irrespective  of  canopy  or  suspension 
line  elasticity,  etc.  However,  the  canopy  and  suspension  line  elasticity,  canopy  mass 
and  drag,  etc.,  all  determine  the  nature  of  the  force  vs.  time  history. 

Examination  of  Fig.  14,  reproduced  from  (22],  clearly  shows  this  effect.  At  a 
given  Froude  number  and  mass  ratio,  the  full-scale  parachute  with  the  4.1  oz/yd^ 
fabric  has  the  largest  non-dimensional  opening  time,  but  the  smallest  non- 
dimensional  opening  shock.  The  V4-scale,  4.1  oz/yd^,  parachute  has  the  shortest 
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rip.  H:  Non-»linH*nsioiial  ojx’niitg  shock  and  non-dimensional  opicnlng  lime  T 
for  three  C-9  lyj'e  [larachtiles,  mass  ratio  0.13. 

(Figure  taken  from  Ref.  22) 


opening  times  over  the  complete  Froude  number  range,  and  tends  to  have  the 
greatest  opening  shock.  The  V4-scale,  0.94  oz/yd^,  parachute  has  intermediate 
values  of  opening  time  and  opening  shock.  Obviously,  canopy  relative  stiffness 
index  affects  this  in  some  manner,  and  produces  the  trends  shown. 

Based  on  the  above,  it  was  reasoned  that  the  product  of  XpX  from  Fig.  14  might 
be  suitable  as  a  correlating  parameter,  where  the  relative  stiffness  index  determines 
the  value  of  x,  as  in  Fig.  8.  The  product  XpX,  called  here  "modified  impulse,"  is 
obviously  related  to  the  total  deceleration  impulse.  The  total  deceleration  impulse 
must  be  some  function  of  the  mass  ratio  and  Froude  number,  since  these  quantities 
are  related  to  the  total  momentum  decrease  the  payload  must  exjjerience. 

To  test  this  hypothesis,  the  product  of  jjeak  opening  shock  and  opening  time, 
XpX  (non-dimensionalized),  was  plotted  as  a  function  of  Froude  number  using  the 
data  from  (9,22).  The  results  are  shown  in  Fig.  15.  Data  for  both  mass  ratios  of  0.13 
and  0.20  are  shown  on  the  same  figure.  The  data  correlation  is  seen  to  be  fair  and 
indicates  that  the  hypothesis  seems  reasonable.  However,  more  work  in  this  area  is 
necessary. 

Two  points  need  to  be  made  regarding  Fig.  15.  First,  all  the  data  points  from 
[9,22]  are  not  on  this  plot  because  it  is  sometimes  difficult  to  tell  from  the  [9,22]  plots 
which  value  of  Xp  goes  with  w’hich  value  of  x  when  several  points  are  plotted  near 
the  same  Froude  number.  However,  Fig.  15  contains  enough  points  over  a 
sufficient  range  of  values  to  define  the  trend.  Second,  there  is  some  difference 
between  the  values  at  mass  ratio  0.13  and  mass  ratio  0.20,  even  though  a  single 
curve  is  faired  through  the  data.  This  probably  reflects  the  theoretical  trend  that 
occurs  if  Figs.  2  and  3  are  combined.  Multiplying  the  data  from  Figs.  2  and  3  together 
gives  the  trends  showm  on  the  theoretical  plots  in  Fig.  15.  Note  that  the  curve  is 
slightly  higher  for  mass  ratio  0.2  than  it  is  for  0.13. 

The  important  fact  here  is  that  the  experimental  trend  for  XpX  vs.  Froude 
number  is  of  the  same  slope  and  approximately  the  same  magnitude  as  predicted  by 
theory. 
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Use  of  Modified  Impulse  Parameter 

The  marmer  in  which  Fig.  15,  the  modified  impulse  graph,  can  be  used  is  as 
follows.  For  a  given  fabric  and  size  of  a  full-scale  or  model  C-9  parachute,  the 
relative  stiffness  index,  is  calculated  using  the  drop  conditions,  dimensions  and 
fabric  properties.  Fig.  8  is  used  to  find  the  parachute  opening  time.  It  should  be 
possible  to  use  this  graph  to  estimate  the  opening  time  of  other  solid  flat  circular 
parachutes  as  well,  as  long  as  the  design  does  not  deviate  too  much  from  the  C-9. 
With  the  Froude  number  and  mass  ratio.  Fig.  15  is  then  used  to  find  the  modified 
impulse.  The  modified  impulse  value  can  then  be  used  to  obtain  the  non- 
dimensional  peak  opening  shock  by  dividing  the  modified  impulse  by  the  non- 
dimensional  opening  time. 

In  this  manner,  the  effect  of  parachute  canopy  stiffness  on  both  opening  time 
and  opening  force  is  taken  into  account.  The  graphs  apply  only  for  mass  ratios  in 
the  vicinity  of  0.13  to  0.20,  however. 

Procedures  for  opening  shock  for  the  T-10  canopy  and  other  canopies  would 
be  similar  to  those  followed  here  for  the  C-9  canopy,  except  that  experimental  data 
over  a  wide  range  of  model  sizes  and  canopy  fabrics  would  be  needed  to  extend  the 
range  for  relative  stiffness  index  shown  in  Fig.  10.  Sufficient  data  on  opening  force 
for  the  T-10  parachute  are  available  from  [27]  to  apply  this  method  over  a  narrow 
range  of  relative  stiffness  indices  for  the  two  canopy  fabrics  used  for  this  design. 
However,  the  time  period  available  for  this  study  did  not  permit  application  of  the 
method  to  the  force  data  for  the  T-10  parachute. 
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OPENING  CHARACTERISTIC8 
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CLOSURE 


Recommendations  for  Future  Research 

The  current  work  has  brought  out  a  number  of  areas  where  further  research 
is  needed.  Any  evaluation  of  experimental  data  is  somewhat  hampered  if  some 
simplified  theories  are  not  available  to  guide  in  interpreting  test  data.  The 
following  quote  taken  from  [30]  aptly  depicts  the  situation  regarding  prediction 
methods: 

The  difficulty  in  predicting  the  performance  of  a  parachute  lies  in  the  complex  interaction 
between  the  porous  canopy  and  the  surrounding  flow  field.  The  parachute  inflation  process 
involves  the  uruteady,  viscous,  compressible  flow  about  a  porous  body  that  undergoes  large 
shape  changes.  Moreover,  this  body  is  composed  of  nonlinear  materials  with  complex  strain, 
strain  rate,  and  hysteresis  properties.  Thus  it  is  not  surprising  that  a  rigorous  analysis  of  the 
Navier-Stokes  equations  for  the  unsteady  flow  about  an  inflating  parachute  presents  a 
formidable  challenge  to  existing  computational  capabilities.  A  recent  study  at  Sandia 
National  Laboratories  concluded  that  an  axisymmetric  flow  field  solution  for  a  typical 
weajwns  parachute  would  require  3300  hours  on  a  CDC  7600  computer  or  330  hours  on  a  CRAY 
system.  Hence,  there  exists  an  urgent  need  for  a  dependable  intermediate  theory  useful  for 
parachute  design. 

The  theories  of  Lingard  and  Fu  referenced  in  the  present  report  seem  to 
provide  a  reasonable  compromise  in  accuracy  without  being  excessively  unwieldy. 
However,  proper  verification  of  these  parachute  opening  theories  would  require  a 
wide  range  of  tests.  The  scatter  in  test  data  creates  this  need.  Most  of  these  tests 
would  have  to  be  conducted  outside  the  practical  range  of  mass  ratios  and  Froude 
numbers. 

Conversely,  at  Froude  numbers  on  the  order  of  four  to  six  where  opening 
shock  loads  are  reasonable,  extrapolarion  of  Fig.  3  or  4  shows  almost  no  effect  of 
mass  ratio  on  non-dimensional  opening  farce  For  this  reason  further  tests  need  to 
be  conducted  at  high  Froude  numbers  over  a  wide  range  of  mass  ratios. 


Use  of  Theories  to  Predict  Canopy  Stiffness  Effects 

It  would  be  most  valuable  to  develop  a  canopy  opening  dynamics  theory  that 
could  include  simulated  canopy  stiffness  effects  in  such  a  way  as  to  confirm  the-  test 
data  trend  shown  in  Fig.  8.  It  appears  that  the  theory  developed  by  Fu  could  be  used 
to  do  this  if  canopy  stiffness  effects  were  simulated  by  using  various  input  values  for 
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initial  semiaperture  angle,  Po  (see  Fig.  1).  Heinrich's  stiffness  index  used  the 
inverted  hang  test  of  the  canopy  and  an  experimental  measurement  of  Dmax  to 
represent  canopy  stiffness.  The  Dmax  and  stiffness  index  values  should  be  reflected 
in  the  initial  semiapertural  angle  P©  some  manner.  Fu's  theory  also  allows  part 
of  the  areal  density  effect  of  the  fabric  to  be  reflected  by  making  an  input  change  to 
canopy  mass.  Some  of  the  effect  of  canopy  stiffness  could  be  accounted  for  by  an 
adjustment  to  suspension  line  elasticity  input,  but  this  does  not  duplicate  the  main 
effect  of  canopy  stiffness. 

Fu  applied  his  theory  only  to  aircraft  drop  test  data  of  full-scale  28  ft  D©  C-9 
parachutes.  The  effect  of  simulated  stiffness  should  be  incorporated  into  the  theory 
and  then  applied  to  Lee's  drop  test  data  for  the  V2-scale  and  V4-scale  C-9  models, 
and  to  the  full-scale  C-9  parachutes  that  were  tested  with  the  heavier  fabrics.  Some 
type  of  correlation  between  a  simulated  initial  semiaperture  angle  and  a  modified 
hang  test  of  the  canopy  should  be  attempted. 


Additional  Verification  of  Relative  Stiffness  Index  and  Modified  Impulse 

Only  for  the  C-9  parachute  is  extensive  test  data  available  for  various  canopy 
sizes  and  fabrics.  The  procedure  used  to  generate  the  opening  time  vs.  relative 
stiffness  index  correlation  plot  (Fig.  8)  should  be  applied  to  another  parachute  tv'pe  as 
well.  Scale  models  of  the  T-10  parachute  would  be  a  likely  candidate  for  this  since  so 
much  full-scale  drop  data  is  already  available  for  this  parachute.  More  work  also 
needs  to  be  done  to  verify  the  modified  impulse  correlation  concept.  The  T-10 
parachute  would  be  a  likely  candidate  for  this  study  al:;o. 


Relative  Stiffness  for  Other  Canopy  Designs 

The  relative  stiffness  index  as  used  in  its  present  form  (Eq.  (2))  has  only  been 
validated  for  the  solid  flat  circular  parachute.  It  would  be  desirable  to  do  research  to 
extend  this  concept  to  more  complicated  geometries  such  as  ribbon,  ringslot,  or  cross 
parachutes. 
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Conclusion 


Theories  for  parachute  opening  dynamics,  dimensional  analyses,  and  test  data 
have  been  reviewed  to  determine  the  effects  of  parachute  canopy  scale  and  stiffness 
on  the  opening  characteristics.  Based  on  this  work,  a  new  relative  stiffness  index  for 
parachute  canopies  has  been  developed. 

It  has  been  shown  that  this  index  gives  good  data  correlation  for  parachute 
opening  times  for  solid  flat  circular  canopies  over  a  wide  range  of  canopy  diameters 
and  fabric  areal  densities.  A  modified  impulse  has  been  developed  based  on  this 
stiffness  index  that  gives  fair  to  good  correlation  for  parachute  opening  shock.  The 
relative  stiffness  index  will  make  it  possible  to  use  scale  model  wind  tunnel  and 
drop  test  data  to  predict  prototype  parachute  performance  with  more  confidence 
than  in  the  past. 

Further  research  needs  to  be  done  to  adapt  the  stiff  ness  index  to  computer  programs  to 
simulate  the  effect  of  canopy  stiffness  in  the  theory.  The  modified  theory  then  needs  to  be 
compared  with  experiments.  Further  work  also  needs  to  be  done  to  adapt  the  relative  stiffness 
concept  to  other  types  of  parachute  canopies. 
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